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INTRODUCTION 


Residual stresses are stresses that remain in a body after all external stresses are removed. 
They may either lower the usable strength, or may be beneficial and raise the strength. Residual 
stresses can arise from manufacturing processes or from non-uniform temperature changes. They 
result from stress levels causing plastic deformation of the material, which is then elastically 
unloaded, but prevented from returning to its initial state by interior or exterior structural restraints. 
The welding process, in particular, creates residual stresses. These stresses have been related to 
cracking in repair welds which has been observed in A1 2195, the lightweight aluminum-lithium 
alloy to be used in the Light Weight Space Shuttle Tank. 

The mechanism that results in residual stresses in the welding process starts with the 
deposition of molten weld metal which heats the immediately adjacent material. After the 
solidification of weld material, normal thermal shrinkage is resisted by the adjacent, cooler 
material. When the thermal strain exceeds the elastic strain corresponding to the yield point stress, 
the stress level is limited by this value, which decreases with increasing temperature. Cooling then 
causes elastic unloading which is restrained by the adjoining material. Permanent plastic strain 
occurs, and tension is caused in the region immediately adjacent to the weld material. Compression 
arises in the metal farther from the weld in order to maintain overall static equilibrium. Subsequent 
repair welds may add to the level of residual stresses. 

The level of residual stress is related to the to onset of fracture during welding. Thus, it is of 
great importance to be able to predict the level of residual stresses remaining after a weld 
procedure, and to determine the factors, such as weld speed, temperature, direction, and number of 
passes, which may affect the magnitude of remaining residual stress. The purpose of this project 
was to develop a simple model which could be used to study residual stress. It was hoped to use 
traditional analytical modeling techniques so that it would be easier to comprehend the effect of 
these variables on the resulting stress. This approach was chosen in place of finite element methods 
so as to facilitate the understanding of the physical processes. The accuracy of the results was to be 
checked with some existing experimental studies giving residual stress levels found from X-ray 
diffraction measurements. 


THEMODEL 


The structure which was modeled was a thin plate of infinite width. It was considered to have 
a weld of uniform thickness deposited along a central line. The plate was assumed to have uniform 
thickness and uniform physical properties, although the yield stress, elastic constants, and 
coefficient of thermal expansion vary with temperature. The coordinate system is shown in Figure 
1. The x-coordinate is along the width of the plate perpendicular from the weld, y is the weld 
direction, and z is the thickness direction. 

The analysis of the mathematical model was composed of a thermal analysis to determine the 
temperature distribution, and a subsequent stress analysis, using the resulting temperature 
distribution. 


THERMAL ANALYSIS 


The heat supplied by a welding arc produces complex thermal cycles in the weld zone. It was 
assumed that the heat source was a point source, and all losses were by conduction into the plate. 
The temperature was assumed to be unchanged at a large distance from the weld. The basic 
equation is: 


dT_ A^T o^T a 2 T\ 

dt (dx 2 dy 2 dz 2 J 
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where 

k = thermal diffusivity = X/cp, X = thermal conductivity, c = specific heat 

It was assumed that the temperature distribution varied only in the x-direction so that the 
equation could be further simplified to 

£T_ / d 2 T 
dt (dx 2 

This equation was solved by a finite difference formulation: 

Iil+1 ~ ^ij - _ -/ ^i+li ~ Tjj _ T ji — Tj-Us 
Axl Ax Ax 


At 


or 


At 


T ij+1 = K ^2 (Ti+1 J " 2Ti i + T i-lj) + T ij 

where the subscript j represents a time 
increment and i represents a spatial nodal 
increment in the x-direction. The model is 
shown in Figure 1. 

Initially, the temperature at the central 
node was set at 600° C, the approximate 
melting temperature of aluminum, and the 
other nodes set to 20° C ambient 
temperature. The above equation was used 
to calculate heat flow over enough time 
intervals until an approximately uniform 
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Figure 1. Finite Difference Model of Plate 
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distribution was reached across the plate. This formulation will approximate the temperature states 
passed through by a welded plate. A total of 27 nodes were used, representing symmetric nodes on 
each side of the center weld. The nodes were spaced at x = weld width = thickness of plate = 1/4” 
A time increment was chosen corresponding to 

a weld speed of v = 8 in/min. Then At = Ay/v 

= l/8"/8 in/min - 1/64 min or approximately 1 
second. However, it was found that a time 
interval of 0. 1 second was required for 

convergence. A value of k = 53mm 2 /sec = .08 

in 2 /sec for 2219 A1 was used. The MATLAB 
program was used to solve the equations, and it 
was found that an approximately uniform 
distribution was reached after 76 time iterations 
or 7.1 seconds when the change in temperature 
at the middle of plate from one time interval to 
the next was less than 5%. The temperature 
distribution for eight time intervals is shown in 
Figure 2. These values were used for the initial 
values in the following stress analysis. 
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Figure 2. Results of Thermal Analysis 
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STRESS ANALYSIS 


Residual stresses were calculated using a model which represents the plate as a series of 
unconnected bars, undergoing only axial strain. Eight time intervals were selected, and each bar 
was assigned a temperature distribution based on the thermal analysis above. The mechanical 
properties were assumed to vary with temperature, and values for yield stress, coefficient of 
thermal expansion, and modulus of elasticity were used for 2219 aluminum alloy based on Ref. 2 
It was assumed that an idealized stress-strain relationship holds, with a constant yield stress 
for each temperature level. Time intervals were taken close together near the onset of welding and 
further apart as the time approached the state of uniform temperature distribution . The effects of 
transverse stress were ignored and uniform longitudinal strain across the plate was assumed. 


The solution procedure was based on one developed in Ref 1 and is as follows: At time 0 the 
thermal stress for the i-th axial element is calculated from Oyj = EjCtjATj, where Ej is the modulus 

of elasticity, cq is the coefficient of thermal expansion, and ATj is the temperature change from the 
previous time interval (or the initial temperature for the first time interval). The axial force over the 
entire width of the plate is then calculated and averaged to produce an average stress in each 
element. In order to maintain static equilibrium, this stress is then subtracted from the thermal 
stress to produce a net stress level. If this value exceeds the yield stress, it is lowered to that value. 
A new axial force is then calculated over the width, and an check for static equilibrium performed. 
If static equilibrium is not satisfied to within a predetermined tolerance, a new average axial stress 
is calculated and the process repeated until static 
equilibrium is attained. Then, the temperature 
distribution for the next time interval is applied, 
and the calculation process begun again, adding 
the stresses resulting from the previous time 
increment, before matching them against the 
yield stress level. At the conclusion of the 
process, the final stress level at the last time 
increment is the residual stress in the plate. 

The results of the analysis are shown in 
Figure 3. 

A maximum value of 48.5 ksi tensile is 
found at the weld centerline, with a 
compressive value of -6.0 ksi seen at the edge 
of the weld. The compressive stress then 
decreases as the distance from the weld 



increases. 

Comparison with Experimental 
Results 

These results were compared with X-ray 
diffraction residual stress measurements 
previously made on 2195 aluminum alloy. 
These results are shown in Figure 4. 

They reveal a maximum tensile stress of 39.7 
ksi at a distance of 0.42 in from the weld 
centerline, and a maximum compressive stress 
of -16.7 ksi at a distance of 1.08 in from the 
centerline. They also show a reduction of 
tensile stress within the weld to almost zero at 
the centerline. The analytical model gives a 
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maximum tensile stress of 48.7 ksi at the centerline and a maximum compressive stress of -6.0 ksi 
at a distance of 0.5 inches from the centerline. These results compare reasonably well in 
magnitude, although the experimental results show stress maxima to be located farther away from 
the centerline. The analytical results show the same shape of stress distribution. 

While this correlation is encouraging for the simple model used, it was hypothesized that an 
improved model, incorporating the effects of transverse stress, might show the decrease in tensile 
stress within the weld seen in the experimental study. It also would provide a value for the 
transverse stress, which may be related to the onset of fracture observed in repair welds. The 
supposition was that as the weld area shrinks, a restraint to transverse contraction is provided by 
the adjacent cooler material. This restraint causes a transverse tensile stress which is associated 
with a longitudinal compressive stress, through the Poisson effect This stress would tend to lower 
the tensile longitudinal stress previously modeled. 

STRESS ANALYSIS INCLUDING TRANSVERSE STRESS 

Each of the elements in Figure 1 is assumed to be in a state of plane stress. It is assumed that 

the longitudinal strain 8yi and the transverse stress is constant across all elements i. The stress 
strain relations in the i-th element are (Ref. 3) 

^xi = 7T"[°x — Vj oFyd + cijTj + 6 x j+ Ae x j 


T r P P 

Ey = - Vi + Oyi] + CtiTj + £yj+ Eyi 

JLsj 

where e^ is the transverse strain, £ y is the constant longitudinal strain, e^ and £ P j are the plastic 

strains which accumulated prior to the current time interval, As^and A£ P j are the incremental 

plastic strains developed in the current time interval, a x is the constant transverse stress, Oy, is the 

constant longitudinal stress, and v; is Poisson's ratio. Static equilibrium across the width of the 
plate is given by the equation 

2a y jAxj = 0 

and the condition restricting the total transverse displacement to zero is given as 

2£ x iAxi = 0 

where the summations are over all the elements, and Axj is the width of the i-th element These 

equations can be written as a series of linear equations in the unknown quantities a y j , e x i e y . 
At the first time interval, t = 0 the right-hand side of these equations contains the known quantities 

ctiTi , where Ti are again taken from the thermal stress analysis. Elastic strains are assumed, so 

that the plastic strain terms c P and As P are zero. The equations are then solved, and if a stress level 
exceeds the yield stress, it is equated to the yield stress. The equations are then solved again, but 
the yield stress for that element is now known, and a plastic strain is now introduced as an 
unknown, replacing the role of the elastic stress component. The solution of this set of equations 
yields the plastic strain increment for the first time interval. At the next time period the previous 

plastic strain increment becomes the e P term and is placed on the on the right hand side of the 
equations. The procedure is repeated for each time interval. The solution at the final time interval 
provides the residual stresses in the plate. The results are shown in Figure 5. 
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Longitudinal and Transverse stress 


They show a similar distribution to the 
previous analysis. The maximum tensile 
stress was reduced about 30% to 34 ksi and 
compares more closely to the experimental 
results. The compressive stress level of -6 
ksi is comparable to the previous result. A 
transverse compressive stress of -4.5 ksi 
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experimental result varying from -4.6 to 
-13.1 ksi located from 1 in to 3 1/2 inches 
from the weld centerline. 
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Figure 5. Results of Transverse Stress Analysis 


Future Directions 


This project did succeed in reproducing qualitative residual stress characteristics similar to 
expected and observed results. The maximum value of tensile residual stress agreed within 12% of 
the experimental values, while the compressive stress levels were approximately 1/2 of the 
experimental results. The dip in tensile values in the interior of the weld was not reflected in the 
model. The following areas for further research in this area should be pursued: 

The thermal analysis can easily be coupled with the stress analysis for greater computational 
efficiency. In this way a solution would flow directly from the assumed temperature distribution to 
the final residual stresses. 

2. Refinement of the model within the weld area should be made to include smaller elements 
and changed properties reflecting the changed metallurgical characteristics of the solidified weld 
material. 


3. A more realistic thermal model reflecting the continuous heat input of a moving line heat 
source would could have a large effect on the residual stresses. 

u fj transv erse stress model should be improved. The plastic strain portion of the solution 
should be analyzed more carefully, possibly using the Prandtl-Reuss flow rule. 

5. Critical stress values should be correlated to fracture of welded plates and to further 
experimental results which should be tied to analytical models. 
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